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ABSTRACT 

p^ ' We present a first systematic comparison of supcrluminous Type la supernovae 

£> \ (SNc la) at late epochs, including previously unpublished photometric and spectro- 

(NJ ■ scopic observations of SN 2007if, SN 2009dc and SNF20080723-012. Photometrically, 

lO ' the objects of our sample show a diverse late-time behaviour, some of them fading 

ON , quite rapidly after a light-curve break at ~ 150-200 d. The latter is likely the result of 

flux redistribution into the infrared, possibly caused by dust formation, rather than a 
^+ | true bolomctric effect. Nebular spectra of supcrluminous SNc la arc characterised by 

f^S . weak or absent [Fe in] emission, pointing at a low ejecta ionisation state as a result 

CO ' of high densities. To constrain the ejecta and 56 Ni masses of superluminous SNe la, 

*~H \ we compare the observed bolometric light curve of SN 2009dc with synthetic model 

^ ■ light curves, focusing on the radioactive tail after ~60d. Models with enough 56 Ni 

• ?-H ' to explain the light-curve peak by radioactive decay, and at the same time sufficient 

^^ , mass to keep the ejecta velocities low, fail to reproduce the observed light-curve tail 

of SN 2009dc because of too much 7-ray trapping. We instead propose a model with 
~ 1 M Q of 56 Ni and ~ 2 M© of ejecta, which may be interpreted as the explosion of a 
Chandrasekhar-mass white dwarf (WD) enshrouded by 0.6-0.7 M Q of C/O-rich mate- 
rial, as it could result from a merger of two massive C/O WDs. This model reproduces 
the late light curve of SN 2009dc well. A flux deficit at peak may be compensated by 
light from the interaction of the ejecta with the surrounding material. 

Key words: supernovae: general - supernovae: individual: SN 2009dc - supernovae: 
individual: SN 2007if - supernovae: individual: SN 2006gz - supernovae: individual: 
SNF20080723-012 - radiative transfer. 



1 INTRODUCTION ll99Sl ; |Perlmutter et al. l ll999h . They are considered amongst 

the most promising tools to distinguish between a Cos- 
Type la supernovae (SNe la) are powerful distance mdi- , . , „ , ,, r r „ , „ 
J v ' mological Constant and other forms of Dark Energy 

cators and have been used to infer the accelerating ex- r^; — , 7-^ — n \ — 71lr,r,i A at 1 ot>t t r^ 5 — T^l 

„ , TT . i-rr n (777^1 U7~, — tt- 52 n (Goobar & Leibundgut 2011). Normal SNc fa (Branch ct al. 

pansion of the Universe Kicss ct al. 1998; Schmidt ct al. 1' on k ,, ,, ., , „ ,, . '. ,, . 

1 I199JI ) are excellently suited tor this purpose owing to their 

remarkable homogeneity in peak luminosity and light-curve 

shape. However, not all SNe la are so 'we ll -behaved'. Several 

Based on observations at ESO Parana!, Pros. 281.D-5043, , , , nTTT: ; — — Tl L _„J rz — ~ ~ ; 71 

n«3 n 0798 H 08=; D 0701 subclasses (e.g. II lhppenko et af.l 19921; ILeibundgut et alj 

f E-mail: tauben@mpa-garching.mpg.de 
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Il993l ; iLi et al.r[2003l ; iHowell et"ail 12009 ; iFolev et al.l |20ld : 
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iGaneshalingam et alj 120121 ') are known to defy normalisa- 
tion through the usu al relations between light-curve w idth 
and peak luminosity (|Phillipslll993l : IPhillips et alJl998l ). one 
of them being the group of superluminous Type la SNe. 

Compared to normal SNe la, superluminous SNe la 
are characterised by a bright light-curve peak, a slow light- 
curve evolution during th e photospheric pha s e and moder- 
ately low eject a velocities (IHowell c t al. 2006; Branch 2006; 
Hicken et alJl20qftlYamanaka et al.ll2009l;IScalzo et alj201(j 



Yuan et al.ll2010l : ISilverman et alj|201ll ; iTaubenberger et alj 
2011T ). Modelling suggests ejecta masses far in excess of 



the Chandrasekhar-mass (Mch) limit of non-rotating white 
dwarfs (WDs) and the production of about 1.5 Mq of 
56 Ni, precluding the interpretation of these events as ther- 
monuclear explosions of Chandrasekhar-mass (Mch) WDs. 
For this reason they are commonly ref erred to as 'super- 
Chandrasekhar' SN e la in the literature (|Howell et al.ll2006l; 
iBranch et alll2006h . 

Models of thermonuclear explosions are severely chal- 
lenged by superluminous SNe la. Proposed explanations 
range fro m the explosions of differentially rot ating massive 
WDs (e.g.. lHowell et al.ll2006t ) to WD mergers (|Hicken et all 
I2007T ) with possible interaction between the actual SN ejecta 
and sur rounding circumstellar medium (CSM) left b y the 
merger (|Scalzo et al.l2010l . l20l3 ; lHachinger et alj2012h . and 
therm onuclear explosions i n the degenerate cores of AGB 
stars IjTaubenberger et alj |2011[ ). However, none of these 
models is without problems. In particular, there is no ev- 
idence of ejecta-CSM interaction in the form of narrow 
emission lines in the photospheric spectra of any of these 
SNe. The problems to explain superluminous SNe la with 
thermonuclear explosions even led to the speculation that 
these objects might instead have a core-collapse origin, with 
the early light curve possib ly powered by magnetar heating 
JTaubenberger et al.lfioilh . 

Great power to discriminate different explosion scenar- 
ios has traditionally been ascribed to nebular spectra. Sam- 
pling also the inner ejecta, differences between a thermonu- 
clear explosion and a core-collapse event should become ev- 
ident in nebular spectra, given the entirely different nucle- 
osynthetic footprint. Moreover, from late-time light curves 
a refined estimate of the 56 Ni mass can be made, comple- 
menting the estimates based on the light-curve peak. 

So far, only few late-time observations of super- 
luminous SN e la h av e been pub li shed. For SNe 2007if 



Scalzo et all [2010: 



( Yamanaka et alj|2009t 



Yuan et al 



l201ll ; ITaubenberger et al 



Tanaka et alj 



20ld ) and 2009dc 



"2010l ; ISilverman et all 
2011) nebular spectra 



shown by lYuan et alj (120101 ); ISilverman et al.1 f|201lh and 



ITaubenberger et alj ( 201ll ). These spectra were dominated 
by [Fell] emiss ion lines, [Fein] features were quite w eak. 
For SN 2006gz l|Hicken et alj|2007h iMaeda et all (|2009T ) pre- 
sented a nebular spectrum that seemed to lack completely 
the prominent [Fe 11] and [Fe ill] emission lines shortward of 
5500 A, normally the hallmark features of SNe la at late 
epochs. At the same time, SN 2006gz was unexpectedly 
dim one year after maximum light, indicating an enhanced 
fading after the photospheric phase. A similar trend, though 
less extreme, was reported for SN 2009dc after ~ 200 d 



i; 



ISilverman et al.ll201ll ; ITaubenberger et al.llioTH ). 

In this work additional late-time data of the superlumi- 
nous SNe la SN 2007if, SN 2009dc and SNF20080723-012 



(|Scalzo et alj |2012j) are presented and analy sed. Comple- 
ment ed by the published data of SN 2006gz JMaeda et alj 
2009), a first systematic comparison of these objects dur- 
ing the nebular phase is performed, revealing significant di- 
versity in the spectral appearance and the luminosity and 
decline rate of the radioactive tail of the light curve. The pa- 
per is organised as follows. In Scction[2]the methods used to 
reduce and calibrate the new data are described. Section [3] 
presents the late-time luminosity evolution, Section [3] the 
properties and peculiarities of the nebular spectra. In Sec- 
tion [5] the findings are discussed, synthetic bolometric light 
curves for a set of models are compared to the observations, 
and an attempt is made to propose a uniform model for at 
least some superluminous SNe la. A summary and conclu- 
sions are given in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

The late-time data presented here were obtained with the 
FORS2 instrument mounted at the 8.2 m Very Large Tele- 
scope (VLT) UT-1, with XShooter at the 8.2 m VLT UT- 
2, and with CAFOS at the 2.2 m Telescope of the Calar 
Alto Observatory. The CAFOS and FORS2 data were re- 
duced following standard procedures within IRAF, includ- 
ing bias subtraction and flatfieldi ng. An optimal, variance- 
weighted extraction |Hornel ll986T ) of the spectra (Table [1} 
was performed, and arc-lamp exposures were used to de- 
termine the dispersion solution. The XShooter observations 
(Table [1} were pipeline-processeqj to create linearised, sky- 
subtracted, wavelength-calibrated 2D spectra (for each of 
the UVB, VIS and NIR channels) out of the curved Echelle 
orders of XShooter. The 2D spectra were then optimally ex- 
tracted. Telluric-feature removal and a rough flux calibration 
were accomplished using observations of spectrophotometric 
(UVB, VIS) or telluric (NIR) standard stars. The flux cal- 
ibration of all optical spectra was checked against the pho- 
tometry and adjusted when necessary. Lacking contempora- 
neous photometry, this was not possible for the XShooter 
NIR spectrum, which was instead scaled to match the VIS 
spectrum in the region of overlap. Since this region is small, 
noisy and affected by the steep transmission edge of the 
dichroic, the obtained calibration is estimated to be no more 
precise than to a factor of ~ 2. 



2.1 Host-galaxy subtraction 

2.1.1 SNF20080723-012 

SNF20080723-012 was located within a spiral host galaxy 
(Fig. [U top). Luckily, it was sufficiently off-set from the 
bright spiral arms that a direct point-spread-function (PSF) 
fitting measurement was feasible under very good seeing con- 
ditions. Such conditions (~ 0.5 arcsec FWHM) were met on 
2009 May 22 with the VLT + FORS2. In the images of that 
night a direct PSF-fitting measurement of the SN magni- 
tudes was performed. At the same time, the PSF-subtracted 



1 http://www.eso.org/obscrving/dfo/quality/XShooter/ 
pipclinc/pipe_gen.html 
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Table 1. Spectra of SNc 2007if, SNF20080723-012 and SN 2009dc. 



UT date 



MJD 



SN 



Epoch" Telescope/Instrument 



Setup 



Exposure time 



2008 Sep 24 

2008 Oct 22 

2009 Apr 26 
2009 May 22 
2009 Jun 19 



54 733.2 
54 761.2 
54 947.3 

54 973.9 

55 001.2 



SN 2007if 
SN 2007if 

SNF 
SNF 
SNF 



2010 May 10 55 326.2 SN 2009dc 



358.2 
384.3 
247.9 
272.7 
298.1 
371.6 



VLT + FORS2 
VLT + FORS2 
VLT + FORS2 
VLT + FORS2 
VLT + FORS2 
VLT + XShooter 



300V 
300V 
300V 
300V 
300V 
UVB,VIS,NIR 



2850 s x 2 
2850 s x 3 
2880 s x 4 
2850 s x 3 
2850 s x 3 
1350 s x 4 



a Phase in rest-fra me days with respe ct to B-band maximum [MJD = 54 348.4 for SN 2007if llScalzo et alj|2010l) . MJD 
SNF20080723-012 l lScalzo et alj|2012h . MJD = 54 946.6 for SN 2009dc llTaubenberger et alJl201lM . 



54 680.9 for 
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Figure 1. V-band image of the SNF20080723-012 field (top 
panel) and /-band image of the SN 2007if field (bottom panel), 
taken with the VLT UT-1 + FORS2 on UT 2009 May 22 and 
2008 September 24, respectively. The fields of view are 4.2 X 2.7 
arcmin 2 , and the regions around the SNe are triply enlarged in 
the insets. North is up and east to the left. The local sequence 
stars around SNF20080723-012 (Table [Alj are indicated. 



2.1.2 SN 2007if 

Background contamination is a more severe problem in the 
late-time observations of SN 2007if. Its dwarf host galaxy at 
redshift z — 0.074 has an almost stellar PSF, and the SN has 
no di scernible offset (Fig.[TJ bottom; see also I Childress et al.l 
120111 ). By the time of our observations, the SN had faded 
below the brightness of the host. Accordingly, PSF-fitting 
magnitudes were dominated by host-galaxy light, and the 
extracted spectra showed clear evidence of an underlying 
continuum. 

The continuum in the spectra was removed by rescaling 
and subtrac ting a spectral templa te of the host galaxy as 
presented bv l Childress et al.l 1)20111 ). In doing so, we assumed 
near zero flux in those spectral regions where emission lines 
from the SN are not expected. 

Since no sufficiently deep imaging templates were avail- 
able in our photometric bands to perform a conventional 
template subtraction, we applied what may be called a nu- 
merical host subtraction. To this aim, we derived synthetic 
host magnitudes in the BVRI bands from the host-galaxy 
spectral template of lChildress et al.l (|2011f ). scaled to match 
the g-band magni tude of the 7 if hos t as reported by the 
same authors and IScalzo et al.l (|2010l ) . We then measured 
the 'SN + host' magnitudes in our VLT images using aper- 
ture photometry with a 2.5 arcsec radius, large enough to 
include virtually all the light fr om the host (and the high-z 
background galaxy reported by ICtiildress et alj|2011f ). Sub- 
traction of the synthetic host fluxes from the 'SN + host' 
fluxes finally yielded SN magnitudes in each bando 

Since the method described above is very sensitive to 
errors related to the measurement or calibration of either 
the 'SN + host' or host magnitudes, it provided question- 
able results in bands where the SN contribution was small 
compared to that of the host. We therefore decided to cali- 
brate only the V^-band SN magnitude in this way, since this 
is the band where the strongest emission lines are located in 
nebular SN la spectra at z = 0.074, and where therefore the 
contrast between the SN and the background is best. Our 
B, R and I magnitudes were instead synthesised from our 
background-subtracted, flux-calibrated (with respect to the 



images were used as pseudo-templates to perform a host- 
galaxy subtraction at earlier epochs when the seeing condi- 
tions were poorer. The possible error introduced by residual 
SN flux in the pseudo-templates is estimated to be small 
(<; 0.1 mag) for the VLT photometry of 2009 April 26 and 
negligible for the earlier Calar Alto photometry. 



2 Note that lYuan et alj i2010ri obtained a s omew hat fainter host 
q-ban d magnitude than IScalzo et al.l d2010l ) and IChildress et alj 
J201lh . which would result in brighter SN magnitudes. The dif- 
ference between the two numbers has been taken into account in 
the error assigned to the host magnitude, and propagated to the 
error of the SN magnitudes. 
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Table 2. 5- and if-corrcctcd photometry of SN 2007if and SNF20080723-012. 



UT date 



MJD 



SN 



Epoch" 



B 



V 



R 



Telescope Seeing 



2008 Sep 24 


54 733.2 


SN 2007if 


358.2 


24.89 ±0.65 c 


24.79 ± 0.41 


25.14 ±0.65 c 


24.36 ± 0.65 c 


VLT 


1.0 


2008 Sep 30 


54 739.8 


SNF 


54.8 


20.72 ±0.16 


20.23 ±0.20 


20.01 ±0.10 




CAHA 


1.6 


2009 Feb 20 


54 882.2 


SNF 


187.3 


23.35 ±0.30 


23.16 ±0.20 


23.17 ±0.31 




CAHA 


1.2 


2009 Apr 26 


54 947.4 


SNF 


248.0 


24.18 ±0.17 


23.92 ±0.14 


24.76 ±0.12 


24.02 ±0.17 


VLT 


0.7 


2009 May 22 


54 973.2 


SNF 


272.0 


24.40 ±0.14 


24.28 ±0.12 


25.25 ±0.22 


24.69 ±0.21 


VLT 


0.5 



a Phase in rest-frame days with respect to S-band maximum [MJD = 
for SNF20080723-012 jScalzo et aI1l2012l) 1. b Stellar FWHM (arcsec) 



54 348.4 for SN 2007if l|Scalzo et alj|2010l) and MJD 
c Magnitudes obtained from an integration of the 



54 680.9 



redshift-corrected and flux-calibrated spectrum with Bessell BRI filters. 



l/-band photometry) spectrum of SN 2007if taken during 
the same night. 



2.2 Photometric calibration; S- and A'-correction 

Our single-epoch photometry of SN 2007if was obtained un- 
der photometric conditions, and the zero points were derived 

r— 1 l — ■ — -| 

from a Stetson (2000) standard field observed on the same 
night. For SNF20080723-012 we had four epochs of photom- 
etry, not all of them obtained in photometric conditions. 
Therefore, a sequence of stars in the SN fie ld (indic a ted in 
Fig. Q] top) was calibrated with respect to IStetsonl (|2000f ) 
standards observed during the photometric nights on 2009 
April 26 and May 22 (Table ETJ. The SN magnitudes in 
all individual nights were then determined relative to this 
sequence of stars. 

To correct for deviations of the instrumental filter 
responses fro m the standard Johnson-Cousins systems 
(jBessellll 19901 ) and the non-negligible redshift of our targets, 
S- and A-corrections (Table |A"2)) were derived from our neb- 
ular spectra and applied to the SN photometry. Since our 
photometric observations of SNF20080723-012 started al- 
ready 55 d after the B-band maximum, but we had no spec- 
tra at that epoch, we instead used a spectrum of SN 2007if 
[actually a comb ination of the ±51 d , ±62 d and ±67 d spec- 
tra presented bv lScalzo et al.l (|2010f )] to calculate the S- and 
A-corrections. The S- and A-corrections for the 187 d pho- 
tometry of SNF20080723-012 were then derived by linear 
interpolation between the 55 d and 248 d values. 

The final, S- and A-corrected magnitudes of SN 2007if 
and SNF20080723-012 are given in Tabled along with their 
associated photometric errors. 



3 LUMINOSITY EVOLUTION 

Absolute BVRI light curves of the candidate super- 
Chandrasekhar SNe 2006gz, 2007if, 2009dc and 
SNF20080723-012 are presented in Fig. H and com- 
pared to those of SNe 2003du and 1991T as prototypes of 
a normal and a luminous SN la, respectively. The absolute 
magnitudes shown in that figure were computed using the 
distance moduli and colour excesses reported in Table [3] 

It can readily be seen from the panels of Fig. [2] that for 
all SNe the late-phase B-, V-, R- and /-band light curves are 
quite similar to one another, so that the important trends 
are preserved when moving on to bolometric light curves. We 



therefore limit ourselves to discussing in detail the pseudo- 
bolometric light curves of all the objects, which are pre- 
sented in Fig. 13] They were obtained by transforming the 
absolute magnitudes into monochromatic luminosities at the 
effective wavelengths of the filters, interpolating the spectral 
energy distribution linearly and integrating over wavelength. 
Zero flux was assumed at the integration boundaries (i.e. the 
blue edge of the B band and the red edge of the / band). 
The resulting BVRI bolometric light curves are expected to 
account for ~ 50 - 75 per cent of the total bolometric flux 
at maximum li ght when the IR contribu tes ~ 10 per cent 
in SN 2009dc (|Taubenberger et alJlJOllh and ~ 10-15 per 
cent in normal SNe la (Mezrag et al. in prep.), and the UV 
contribution blueward of th e U band should have dropped 
below the 10 per-cent level (|Silverman et al.l 1201 ll ). At late 
phases, we may miss a significant amount of IR flux, but 
lack of late-time IR photometry of superluminous SNe la 
prevents a more quantitative assessment. 

Focusing on the early light curves (inset in Fig. |3J), 
SN 2003du has the narrowest and least luminous peak 
(Lgyftj ~ 0.8 x 10 43 ergs -1 ) in our sample of objects, in 
line with the expectations for a normal SN la. SN 1991T 
is clearly more luminous and more slowly declining than 
SN 2003du, and interestingly quite similar to the proposed 
super-C handrasekhar SNe 2006gz and SNF20080723-012 
(see also lScalzo et al.ll2012r ). These three SNe reach peak lu- 
minosities of Z/svh/ ~ 1.2 x 10 43 ergs -1 . Finally, SNe 2007if 
and 2009dc outshine all the other SNe, being more than 
twice as luminous as SN 2003du at peak {L^mi ~ 1.9 x 10 43 
ergs ). The early-time pseudo-bolometric light curves of 
these two objects are essentially identical. 

During the early nebular phase, between ~ 100 and 
200 d after maximum, all SNe evolve as expected from their 
behaviour around peak brightness. SNe 2007if and 2009dc, 
which feature the brightest peak, are also most luminous at 
these phases. SNF20080723-012 and SN 1991T lie between 
normal SNe la and SNe 2007if and 2009dc. The light-curve 
decline of all these SNe is quite similar, a bit faster in nor- 
mal SNe la and a bit slower in the superluminous SNe, but 
always somewhat faster than the decay rate of 56 Co, which 
is expected in the case of increasing 7-ray losses. No data 
are available for SN 2006gz during the early nebular phase. 

The simple hierarchical picture that seems to emerge 
from studying the luminosity evolution out to ~ 200 d, how- 
ever, does not easily fit the data thereafter. SNe 2003du, 
1991T, 2007if and SNF20080723-012 still decline with 
a nearly constant slope between 200 and 400 d after 
maximum. Among these SNe the luminosity differences 
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Table 3. Basic properties of proposed supcr-Chandrasckhar SNe la. 



SN 2006gz 



SN 2007if SNF20080723-012 SN 2009dc 



Redshift z^el 

Distance modulus fj, (mag) 

Colour excess E(B — V) (mag) 

Am 15 (B) 



0.0237 a 
35.03±0.06 d 
0.18±0.05 a 
0.69±0.04 a 



0.0742 6 
37.45±0.05 b 

0.079 d > e 
0.71±0.06 d 



0.0745 6 
37.46±0.05 b 

0.064 e 
0.93±0.04^ 



0.0214 c 
34.86±0.08 c 
0.17±0.07 c 
0.71±0.03 c 



a Hick cn et al.l d2007r)- b Helioce ntric redshift from Ha e mission in the host galaxy, d istance modulus fro m zwi . c ISilverman et al,l 
l|201lT ); lTaubenberger et al.l l|201ll l. d IScaIzo etHl l|20ld l. e Galactic colour excess l|Schlegel et al.lll998h . f IScalzo et alj J2012h . 
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Figure 2. Absolut e BVRI li ght curves of the proposed super-Chandrasekhar SNe 2006gz JHicken et al.ll2007l; iMaeda et al.ll2009h , 2007if 
llScalzo et al.ll2010ri . 20 09dc ^Silverman et aIll201ll;lTaubenbereer et al.ll2(5Tlh and SNF2008 q723-012 llScalzo et alj|2012h . the 'classical' 
luminous SN la 1991T l lLira et al.lll998ljAltavilla et al.ll2004h and the normal SN la 2003du llStanishev et al.ll2007h . Large symbols with 
error bars are data newly presented in this work (Table [2}- The phase is given in rest-frame days with respect to B-band maximum. 



observed at earlier times are approximately preserved, 
though SNF20080723-012 now appears a little fainter than 
SN 1991T. SN 2006gz, on the other hand, fades dramati- 
cally s ome time between the peak phase and the nebular 
phase (jMaeda et al.ll2009r ). and a year after the explosion 
it is a factor of ~ 4 less luminous than the normal SN la 
2003du and almost an order of magnitude less luminous 
than SN 1991T. SN 2009dc seems to share the destiny of 
SN 2006gz, though to a lesser extent. Around 200 d after 
maximum its light curve starts to decline more rapidly than 



before (|Silverman et all l201ll ; iTaubenberger et all 120111 ). 
and, while not fading as rapidly as SN 2006gz, at 400 d after 
maximum it is no longer more luminous than SN 2003du. 



4 SPECTROSCOPIC COMPARISON 

Nebular spectra of our sample of superluminous SNe la, 
SN 1991T and SN 2003du are shown in Fig. g] In general, 
the spectral evolution during the nebular phase is slow, but 
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Figure 4. Nebular spectra of SNe 2007if l|Blondin et alj|2012l ; this work), 2009dc JTaubenberger et aljlioill ; ISilverman et alj|201ll: this 
work ) and SNF20080723-012 (this work; cf. Table [J , comp ared to those of the proposed su per-Chandrasekhar SN 2006gz JMaeda et alj 
120091) , the luminous SN la 1991T llGomez fc Loped Il998f) and the normal SN la 2003du llStanishev et al.ll2007h . Phases are given in 
rest-frame days with respect to B-band maximum. Left panel: Spectra ~ 120-160 d after maximum. For presentation purposes, the 
spectrum of SNe 2007if and the + 162 d spectrum of SN 2009dc have been smoothed by 500 and 2500 km s~ 1 , respectively. Right panel: 
Spectra ~ 250-400 d after maximum. The two spectra of SN 2009dc and the spectra of SN 2007if and SNF20080723-012 have been 
smoothed by 300, 1300, 1000 and 1000 kms -1 , respectively. The gap in the +372 d spectrum of SN 2009dc is the region between the 
XShooter UVB and VIS channels, with very low signal in both. 



not negligible. Since the spectra in Fig. [4] have been taken 
at very different epochs, evolutionary effects and intrinsic 
differences between SNe have to be disentangled. This is 
straightforward for those objects for which multi-epoch spec- 
troscopy is available. 

Concentrating first on the early nebular phase (left 
panel of Fig. [4}, we see that the superluminous SNe la 2007if 
and 2009dc are relatively similar to each other and also 
to the normal SN la 2003du. The same [Fein], [Fen] and 
[Co in] emission lines are detected in all three SNe (see e.g. 



iMaeda et aT]|2010l ; iMazzali et"aLll201ll . for a detailed iden- 
tification of emission lines in nebular SN la spectra). Dif- 
ferences are evident in the width of spectral features, with 
many lines being resolved into double or multiple peaks in 
SN 2009dc that are just blended into a single broad fea- 
ture in SN 2003du. Also the ionisation state is noticeably 
different, with weaker [Fein] and more prominent [Fen] 
lines in SN 2009dc compared to SN 2003du. In both re- 
spects, SN 2007if takes an intermediate position between 
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same SNe as in Fig. [2] The peak phase is enlarged in the inset. 
Epochs are given in rest-frame days after B-band maximum. 



SNe 2003du and 2009dc, but in terms of ionisation closer to 
SN 2009dc. 

During the later nebular phase (right panel of Fig. [4j 
the observed differences become more pronounced. In nor- 
mal SNe la such as SN 2003du, the most important evo- 
lution between four months and one year after the ex- 
plosion is the fading of Co emission lines due to the ra- 
dioactive decay - of 56 Co into 56 Fe with i 1/2 ( 56 Co) = 77.2 d 
(jKuchner et al.lll994T ). The same is also observed in the su- 
pcrluminous SNe 2007if and 2009dc, but here a lower ejecta 
ionisation state is seen as an additional effect. The hall- 
mark feature of neb ular SN la spectra, the prominent [Fe ill] 
blend at ~ 4700 A (|Axelrodl Il980l ; ISpyromilio et al.l 1 19921 ; 
iMazzali et al.lll998l ). is merely a stump in SNe 2007if and 
2009dc. Instead, most of the emission blueward of 5500 A 
probably originates from [Fen] transitions. 

In normal SNe la the features between 7000 and 7500 A 
are attributed to forbidden transitions of iron-group el- 
ements (IGEs), notably [Fen] A7155 and [Nin] A7378, 
which often fo rm a double-peaked structure as obse rved 
in SN 2003du (|Maeda et al-lboiot Ffanaka et al.ll201ll ). In 
the nebular spectrum of SN 1991T, otherwise very simi- 
lar to that of SN 2003du, there is only a single, rounded, 
broad peak from 7000 to 7500 A, which may be the ef- 
fect of a d istribution of IGEs out to higher velocities in 
SN 1991T |Mazzali et al.l 1 19951 ). In SNe 2007if and 2009dc 
the emission in that region is stronger than in the previ- 
ously mentioned objects. The SN 2007if spectrum is too 
noisy to study the line profiles in detail, but in SN 2009dc 
there are two pronounced, sharp emission peaks. The blue 
peak can be explained by [Fen] A7155, whereas the red- 
der at ~ 7305 A is not at the right position for [Ni n] A7378 
emission. We instead suggest a significant contribution of 
[Can] AA7291, 7324. These lines are typically weak or ab- 
sent in normal SNe la, but have been identified in pecu- 
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Figure 5. Nebular N IR spectrum of SN 2009 dc, compared to the 
normal SN la 1998bu jSpvromilio et al,l2004l) . Phases are given in 
rest-frame days with respect to .B-band maximum. The spectrum 
of SN 2009dc has been smoothed by 2000 kms , and regions 
with strong telluric features are omitted. 



liar SNe la, in particular in subluminous, 91bg- like events 
(JMazzali et al.lll997l : llvlazzali fe Hachingeilfeoi^ . The fact 
that in SN 2009dc the 7000 to 7500 A region shows an 
unusual triple-peaked structure at phases around day 150 
(Fig. 21 left panel) suggests that [Can] emission might be 
present already at those epochs and strengthen with time. 

The nebular spectrum of SN 2006gz shares similari- 
ties with those of SNe 2007if and 2009dc, especially in the 
strength of the 7000 to 7500 A emission. However, the flux 
in the blue part of the spectrum is strongly suppressed, and 
individual features cannot be identified as a consequence of 
the low signal-to- noise ratio. Note that the SN 20 06gz spec- 
trum as presented here and in lMaeda et al.l (|2009j) is binned 
over 16px/22A. The SNF20080723-012 spectrum, on the 
other hand, is very similar to the spectrum of SN 1991T, in- 
cluding the broad single-peaked emission between 7000 and 
7500 A. It thus deviates strongly from the spectra of the 
other proposed super-Chandrasekhar SNe la. Taking into 
account the light-curve properties discussed in the previous 
section, one may therefore speculate that SNF20080723-012 
might be a 91T-like object rather than a true member of the 
cla ss of 'super-Chandr asekhar' SNe la, as already discussed 
bv lScalzo et~aH (|2012t ). 

At NIR wavelengths (Fig. [5| only a few features can 
be safely identified in the spectrum of SN 2009dc as a 
consequence of the low signal-to-noise ratio. The detected 
[Fen] 1.257, 1.279 /im, [Fen] 1.533 /im/ [Con] 1.547 /im and 
[Fen] 1.644 /im emission lines are characteristic of nebular 
NIR spectra o f SNe la, and promine nt also in the normal 
SN la 1998bu (|Spyromilio et al.ll2004l ) which is included in 
Fig. [5] for comparison. A contribution of [ Sil] 1.646 /.mi to 
the 1.64 /xm line (e.g. IMazzali et al.l 120111 ) cannot be ex- 
cluded in SN 2009dc, though the feature is not stronger 
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Figure 6. Relation between Ami5(B) and the flux ratio of the 
[Fe III] / [Fe II] emission lines near 4700 and 5200 A, determined in 
the nebular spectra closest to day 300 after maximum. The line 
fluxes have been measured from zero intensity through a multiplc- 
Gaussian deblendin g of the 4500—5500 A region as described by 
iMazzali et al.l 1119981 ). Red symbols stand for subluminous SNe la, 
blue symbols for proposed super-Chandrasekhar events, and black 
symbols for normal or 91T-like SNe la. 



than in SN 1998bu wh ere it was explained by [Fell] alone 
(|Spvromilio et al.l 12004 1. An additional emission feature is 
identified in SN 2009dc at ~1.17/im, which has no dis- 
cernible counterpart in SN 1998bu. It is located in a re- 
gion devoid of strong telluric absorptions, and is similarly 
pronounced as other lines discussed before. Hence, we tend 
to believe that it is not a mere reduction artefact, but can- 
not provide a conclusive identification without detailed mod- 
elling. 



5 DISCUSSION 

5.1 Ionisation state in the nebular phase 

The ionisation state of SN la ejecta during the nebular 
phase can be assessed by studying the flux ratio of the 
emission blends close to 4700 and 5200 A. As discussed 
in detail e.g. by IMazzali et al.l (|201lh . the 4700 A blend is 
clearly dominated by [Fe ill] emission, whereas the 5200 A 
blend has a significant contribution of [Fell] next to [Fein], 
making the ratio a sensitive ionisation indicator. In or- 
der to reproduce the large observed A4700/A5200 ratio in 
the normal o r slightly sublumino us SN 2003hv with syn- 
thetic spectra. IMazzali et alj (|2011f ) had to assume a reduced 
central density compared to that of the o ne-dimensional 
Chandrasekhar-mass explosion model W7 (|Nomoto et alj 
ll984T ). leading these authors to speculate about a sub- 
Chandrasekhar-mass or merger origin of SN 2003hv. In an 
analogous way, the low ionisation found in SNe 2007if and 
2009dc might be indicative of high central ejecta densi- 
ties, leading to enhanced recombination. This should be 



a direct consequence of the low ejecta expansion veloc- 
ities determined f r om both early- and late-time spectra 
dScalzo et aill2010l; lYamanaka et al.l 120091 ; ISilverman et al.l 
l201ll ; iTaubenberger et al.ll201ll ; see also Section [4}. 

In Fig. [6] the flux ratio of the emission lines near 4700 
and 5200 A is investigated as a function of Ami5(B) [which 
in normal S Ne la may be t aken as a proxy for the peak lu- 
minosity; see lPhillipsI (|l993l )]. Normal and 91T-like SNe (in- 
cluding SNF20080723-012) show flux ratios between 1.3 and 
1.9, SNe 2007if and 2009dc between 1.0 and 1.1. The only 
other SNe with similarly low flux ratios (and hence ionisa- 
tion) are subluminous 91bg-like SNe. Their nebular spectra, 
however, show a complex structure with very narrow [Fe ill] 
lines superimposed on broad [Fe 11] emission (jMazzali et alj 
1 19971 ; IMazzali fc Hachingerll2012r i. 

The inferred low ionisation is also consistent with the 
likely detection of [Can] AA7291, 7324. The first and second 
ionisation potentials of Ca are lower than those of Fe. Ac- 
cordingly, in regions with significant Fein content, Ca should 
be doubly ionised to almost 100 per cent, which is proba- 
bly why [Call] lines are absent from the nebular spectra of 
normal SNe la. The lower ionisation in superluminous and 
91bg-like SNe la, however, should favour Call as the domi- 
nant ionisation state. Indeed, [Ca 11] emission is observed in 
exactly those objects. 



5.2 Late-time dust formation? 

In Section [3] we discussed the diverse light-curve decline of 
superluminous SNe la at late phases. Once on the radioac- 
tive tail ~50d after maximum light, SN 2007if fades at a 
nearly constant rate for more than 300 d. SN 2009dc shows 
a steeper decline after ~200d, and is about a factor 3 less 
luminous than SN 2007if after one year, though it had the 
same peak luminosity (Fig. [3} . SN 2006gz fades even more 
rapidly, being about a factor 15 less luminous than SN 2007if 
one year after maximum, though the difference at peak was 
merely a factor 1.6. 

Following the reasoning of ITaubenberger et alj (|2011T ) . 
the accelerated decline observed in the BVHJ-bolometric 
light curve of SN 2009dc (and probably also SN 2006gz) 
at those late phases is unlikely to be a true bolometric ef- 
fect. As long as the light-curve tail is powered by radioac- 
tive decay and positrons are fully trapped, the decline may 
slow down (when a longer-lived radioactive isotope starts 
to dominate the energy deposi tion), but not accelerate [cf. 
iRuiz-Lapuente fc Spruitl (|l998T ) for the case of incomplete 
positron trapping]. Instead, the observed luminosity drop 
is probably the outcome of flux redistribution into regimes 
that are not observed, most likely the IR. This could be ac- 
complished by an early IR catastrophe or by dust formation. 

Th ough never ob s erved in SNe la to da te, an IR catas- 
trophe (|Axelrod 1980; Fr ansson et al.ll 19961 ) is an inevitable 
consequence of the decreasing energy deposition by radioac- 
tive decay and the expansion of the SN ejecta. Below a crit- 
ical temperature T c the upper levels of forbidden transitions 
in the optical and NIR can no longer be populated. The cool- 
ing is henceforth dominated by ground-state fine-structure 
transitions of Fe in the mid- and far-IR, accompanied by a 
rapid temperature decrease of the ejecta. An IR catastrophe 
is favoured by low densities (as a consequence of the density 
dependence of T c ). In the high-density ejecta of SN 2009dc 
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Figure 7. Comparison of the late-time spectrum of SN 2006gz with that of SN 2009dc, artificially reddened by E(B — V) = 0.2 mag 
(upper left panel), and that of SN 2007if, artificially reddened by E(B — V) = 0.5 mag (upper right panel; Ry = 3.1 assumed in both 
cases). The residuals after subtracting the spectra (lower panels) hardly show any differences exceeding the noise level. 



an early onset of the IR catastrophe is therefore not ex- 
pected. Moreover, even one year after maximum, when the 
-BVRJ-bolometric luminosity of SN 2009dc is already signif- 
icantly below that of SN 2007if, the spectra are still similar, 
showing prominent [Fen] emission throughout the optical 
and NIR wavelength range. 

Dust formation, on the other hand, is usually associated 
with core-collapse SNe and has rarely been discusse d in the 
context of thermonuclear SNe IjNozawa et al.l 12011! ). In su- 
perluminous SNe la, however, it seems to be compatible with 
observations. Early-time spectra of SNe 2006gz and 2009dc 
show C n lines with a strength unprecedented in thermonu- 
clear SNe. This suggests that at least a moderate amoun t of 
carbon is present in the ejecta [Hachinger et al.l (|20T2r ) re- 
produced the spectral time series of SN 2009dc with carbon 
mass fractions between one and ten per cent]. The carbon 
may give rise to the formation of graphite dust, but can also 
pave the way for dust formation in general through prior CO 
molecule formation, which opens an efficient cooling chan- 
nel, the emission in molecular bands. Dust formation should 
also be promoted by the comparatively high densities that 
the slowly ex panding ejecta of su perluminous SNe retain at 
late phases (jNozawa et al.l 120111 ) . The luminosity drop in 
SN 2009dc comes along with an evolution towards redder 
colours (Fig. [2]), again consistent with dust formation. In- 
terestingly, even the seemingly unique nebular spectrum of 
SN 2006gz with its apparent lack of features in the blue could 
find a natural explanation within this scenario. This is shown 
in Fig. [7J where the late-time spectra of SNe 2007if and 
2009dc have been artificially reddened to simulate the effect 
of more pronounced dust formation in SN 2006gz. Nothing 



can be said about individual features blueward of ~ 6000 A 
since the spectrum of SN 2006gz is too noisy in that region, 
but from the overall spectral shape it is plausible that the 
late-time spectra of SNe 2006gz, 2007if and 2009dc are all 
intrinsically similar, and differ just by the amount of red- 
dening caused by newly formed dust in the ejecta. 

Note that changes in the profiles of nebular emission 
lines, usually another signature of dust formation within 
the ejecta, are neither observed nor expected in SN 2009dc. 
If dust formation occurs, it would be most effective in the 
carbon- and silicon-rich zones, which, according to the ve- 
locity ev olution of the C H A6580 an d Sill A6355 lines pre- 
sented bv lTaubenberger et al.l (|201 lh , should be located out- 
side ~ 6000 kms -1 . The [Fell] emission in the nebular spec- 
trum of SN 2009dc, however, comes from below 6000 kms" 1 ; 
emission above that velocity would result in the nebular lines 
being too broad. Accordingly, all parts of the nebular [Fe n] 
emission would be attenuated by the same amount. 

Under the premise that dust formed in the ejecta of 
SN 2009dc, some of its properties can be estimated. To cal- 
culate the dust mass resulting in an extinction Ay = 0.93 
mag [E(B — V) — 0.3 mag, Rv — 3.1] on day 372 after peak, 
we have assumed that the dust is located in an (infinitesi- 
mally) thin shell at radius R. For this geometry, the optical 
depth that corresponds to this value of Ay (i.e. ry = 0.86) 
can be expressed as 



ry 



)v(a)ira 2 M dust 
4tyR 2 m Erain 



(1) 



where a is the dust-grain radius, m gra i n the dust-grain mass, 
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Figure 8. Nebular (+372 d) spectrum of SN 2009dc, corrected for the Galactic and host-galaxy reddening of E(B — V) = 0.17 mag. 
The solid lines are blackbody spectra for temperatures Tbb between 400 and 800 K, assuming emission from a thin dust shell at 8000 
kms -1 . The lighter shaded regions show the effect of varying the shell location between 6000 and 10000 kms - 1 . It appears that a dust 
temperature of > 700 K is excluded by a comparison to the observed NIR spectrum of SN 2009dc. 



Mdust the total dust mass, and Qv(a) the extinction effi- 
ciency in the V band. It follows that 



w !6tT 2 Pgrain O 

Mdust = -5- it „ , r T V , 
3 Qv(a) 



(2) 



very similar to eq. 6 of lLucv et al.l (|1989f ). but with a differ- 
ent prefactor reflecting the different geometry (a thin shell 
here versus a homogeneous sphere there) . For small grain s 
Qvip) is approximately proportional to a (|Lucv et al J 19891 1. 
and hence the expression for Mdust becomes largely inde- 
pendent of a. Assuming a dust-grain density of p gra i n = 2.9 
gem - ' 5 (appropriate for a mixture of graphite and silicateqj), 
a grain s ize of 0.01 /an, and adopting a Qv of 0.07 from 
fig. 4a of iDraine fc Led l|l984f l. we have evaluated Eq. [2] for 
different radii of the dust shell, between 6000 and 10 000 
kms -1 in velocity space. The resulting dust mass ranges 
from ~lx 10~ 4 to ~ 4 x 10 -4 Mq, similar to what has been 
reported for s ome core-collapse SNe in the literature (e.g. 
for SN 2004et; iKotak et al.ll2009h . 

The newly formed dust should also manifest in the 
emission of a thermal continuum reflecting the temperature 
of the dust grains. The emission would likely peak in the 
mid IR, but depending on the dust temperature the NIR 
regime may also be affected by the Wien tail of the black- 



3 When conditions are appropriate for dust formation in superlu- 
minous SNe la, the presence of carbon and also silicon may allow 
for the formation of carbonaceous dust as well as silicates (typical 
densi ties 2.2 and 3.5 gem -3 , respectively; IWeingartner fc Drainel 
I2001T) . Given our ignorance of the exact dust composition we as- 
sume a 1:1 mixture in our calculation; the resulting uncertainty 
in the estimated dust mass is less than ±25 per cent. 



body spectrum. To derive limits on the dust temperature in 
SN 2009dc, we have calculated blackbody curves for differ- 
ent temperatures Tbb and radii R of 6000 to 10 000 kms -1 
in velocity space. The luminosity can be expressed as 



in 2 R 2 B v (T B 



(3) 



In Fig. [5] the derived blackbody curves are overplotted on 
the optical-through-NIR spectrum of SN 2009dc taken 372 d 
after maximum light. The figure suggests that dust temper- 
atures of 700 K or more are excluded by the continuum level 
of the NIR spectrum, whereas ^ 600 K may be acceptable. 
Such temperatures correspond to blackbody spectra peak- 
ing in the mid IR (at > 4 /mi). However, the uncertainties 
in this estimate are large: not only the radius of the dust 
shell is poorly constrained, also the flux calibration of the 
NIR spectrum of SN 2009dc is uncertain by a factor of ~ 2 
(see Section [2]). 

In the end, mid-IR observations of superluminous 
SNe la during the nebular phase will be the only way to 
prove or disprove dust formation on solid grounds. If dust 
forms, excess emission in the mid IR with a thermal spectral 
energy distribution should be detected. 



5.3 Bolometric light-curve models 

The bolometric light curves presented in Section[3]can serve 
as a benchmark for models proposed for superluminous 
SNe la. To enable such a model-data comparison, we have 
computed synthetic bolometric light curves for a number of 
mod els using the Monte- Carlo radiative-transfer code AR- 
TIS ()Kromer fc Simll2009l ). Full spectral calculations at late 
phases are not possible with ARTIS, since non-thermal pro- 
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Figure 9. Synthetic bolomctric light curves (black solid lines) 
for different models compared to th e observed UBVRIJHK- 
bolometric light curve of SN 2009dc l lTaubenberger et alj|201ll . 
blue data points, assuming a B-band rise time of 23 d). The con- 
tributions of 7-rays (red solid lines) and positrons (purple solid 
lines) to the synthetic bolomctric light curves are shown individ- 
ually. 

Top panel: The 09dc-cxp model of lHachinger et alj d2012r). 
Middle panel: The 09dc-int model of lHachinger et al.l J2012T) . 
Bottom panel: 09dc-tail (~ 2 M Q of ejecta, ~ 1 M Q of 56 Ni). 



cesses and a detailed treatment of collisions and forbidden 
transitions are presently not implemented. However, for the 
late-phase bolometric light curve the deposition of 7-rays 
and positrons is the only relevant physical process. This 
is treated with sufficient accuracy in ARTIS. We have run 
our c alculations with a simple grey optical opacity (|Siml 
[2007]), which is a good approximation for the late bolomet- 
ric light curve when ti me-dependent effects are negligible 
l|Cappellaro et al.lll997l ). The rise time, the detailed shape 
of the light-curve peak and the peak luminosity, however, 
may be altered by this approximation. 

For our li ght-curve calcula t ions w e have used the same 
models as in lHachinger et al.l (|2012f ). These authors used 
early-time spectra of SN 2009dc to perform a tomographic 



study of the outer layers of the ejecta for three explosion sce- 
narios. The first type of model investigated is the detonation 
of a rotating supermassive WD with 2 Mb, represented by 
an AWD3det-based density structure IjPfannes et al.ll2010l . 
Fink et al. in prep.). The two alternative scenarios are a 
possible core-collapse origin of SN 2009dc, represented by 
an empirically derived exponential density structure with 3 
Mq ('09dc-exp'), and a model where part of the luminos- 
ity is assu med to be generated by CSM interaction, based 
on the W7 |Nomoto et al.lll984l ) density structure with 1.4 
Mq of ejecta but a rescale d kinetic energy ('09dc- int'). For 
the AWD3det-based model. lHachinger et al.l (|2012l ) were not 
able to obtain a consistent solution: the line velocities in the 
synthetic spectra were too high, and the inferred abundance 
structure was inconsistent with the kinetic energy of the 
underlying density profile. Accordingly, this model was dis- 
favoured on the basis of the early spectral evolution. For 
the 09dc-exp and 09dc-int models, on the other hand, more 
plausible solutions could be obtained. 

In their tomographic study, lHachinger et al.l (J2012J) 
used the luminosity mostly as a free parameter. Synthetic 
bolometric light curves thus provide an additional constraint 
on the structure of the SN by constraining the total energy 
of the exp losion. To compute bol ometric light curves for the 
models of lHachinger et al.l (|2012j ). an assumption had to be 
made for the inner core of the ejecta (below ~ 4000 km s _1 ) 
whose composition was not constrained by the tomography. 
Here we kept the mass fraction of stable iron the same as 
in the innermost shell studied by Hachinger et al., and filled 
the rest with 56 Ni. 

Fig. [9] shows synthetic bolometric light curves of 
the 09dc-exp model (top panel) and the 09dc-int model 
(middle panel), compared to the UBVRIJHK pseudo- 
bo lometric luminosity evoluti on of SN 2009dc as presented 
by iTaubenberger et al.l (|201ll ). Clearly, the 09dc-exp model 
manages to reproduce the peak luminosity of SN 2009dc 
(which is a consequence of the large 56 Ni mass of almost 
1.6 M0), but it is also evident that the radioactive tail is 
much too bright in the model owing to the large total mass 
of 3 Mq that leads to strong 7-ray trapping. This is fur- 
ther illustrated by the red and purple lines in Fig. |5J which 
give the individual contributions of energy deposited by 7- 
rays and positrons, respectively. The intersection thus marks 
the epoch where the 7-ray trapping has decreased to ~ 3 
per cent, which is approximately the fraction of 56 Co de- 
cay energy carried by positrons. In 09dc-exp this level is 
reached very late, ~440d after the explosion. In contrast, 
the 09dc-int model fails to reach the observed peak lumi- 
nosity of SN 2009dc by quite some margin, though in this 
case by purpose: the remainder is supposed to be gener- 
ated by ejecta-CSM interaction, and not included in the 
synthetic bolometric light curve that shows only the con- 
tribution from radioactive decay. However, also the light- 
curve tail of 09dc-int seems much too dim at phases between 
100 and 200 d when the spectra with their lack of pseudo- 
continuum emission do not support si gnificant ejecta-CSM 
interaction (|Taubenberger et al.luOllf ). The model, as con- 
sidered here, has only Chandrasekhar mass, i.e. it repre- 
sents the exploded W D without any swept-up material (cf. 
lHachinger et al.ll2012T ). The low mass leads to inefficient 7- 
ray trapping (reaching 3 per cent already 280 d after the 
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10 15 

v (10 3 km s" 1 ) 

Figure 10. Density profiles of the models shown in Fig. [9] eval- 
uated at a reference time of 100 s after the explosion. 



explosion), insufficient to explain the observed luminosity 
between 100 and 200 d. 

We now test whether the shortcomings of 09dc-int in 
fitting the light-curve tail before the break at ~ 200 d can 
be overcome by adding more mass. To this end, we have 
constructed a new model, '09dc-tail' (Table [4] and bottom 
panel of Fig.[§]). Its density profile p(v) is a moderately steep 
exponential for v ^ 2500 kms -1 , but flattens below that ve- 
locity. This results in a total ejecta mass of ~ 2 Mq , and a 
kinetic energy of 1.2 x 10 51 erg. The density profile of 09dc- 
tail is shown in Fig. 1101 together with those of 09dc-exp 
and 09dc-int. The 56 Ni mass of 09dc-tail was adjusted to be 
1 Mq, more than in 09dc-int (~0.6 Mq) but significantly 
less than in 09dc-exp (~ 1.6 Mq). Thanks to the fairly steep 
density profile and the high ejecta mass the 7-ray opacity 
is large, resulting in 3 per cent 7-ray trapping at 390 d af- 
ter the explosion. The bolometric light curve of the model 
provides a very good fit to the data of SN 2009dc between 
60 and 190 d after the explosion. The model remains more 
luminous than SN 2009dc at later epochs, but as discussed 
in the previous section this may be due to a significant IR 
flux that was missed when the pseudo-bolometric light curve 
of SN 2009dc was constructed from the observed bands. As 
09dc-int, the model does not reach the observed peak lumi- 
nosity of SN 2009dc. It peaks at a flux 30-40 per cent too 
low. 



5.4 Putting the pieces of the puzzle together 

In the previous sect i on we have seen that the 09dc-exp model 
of lrlachinger et al.l (|2012r ) significantly overestimates the ra- 
dioactive tail with respect to the observations of SN 2009dc. 
Though this is just a single model, the consequences of this 
finding are far-reaching. This is because 09dc-exp is quite 
generic for models that try to explain the peak luminosity 
of superluminous SNe la by radioactivity (and hence need 
> 1.5 M© of 56 Ni), and at the same time keep the kinetic 
energy per mass low [which requires a significant amount 
of unbu rned material and a t otal e jecta mass of > 2.5 Mq; 
see also iTaubenberger et al.l (J201l|)]. In these models there 
is no obvious way to avoid the strong 7-ray trapping seen in 



09dc-exp. The consequence is that probably all models with 
excessively large 56 Ni masses will fail. 

If this is true and the luminosity of superluminous 
SNe la around peak does not solely come from radioactive 
decay, other processes have to contribute. Especially ejecta- 
CSM interaction has to be considered as a possibility, despite 
the known caveats that the CSM would have to be H- and 
He-free in order to avoid narrow lines of those elements in 
the spectra of superluminous SNe la, and that even with a 
C/O CSM it is not evident that telltale spe ctroscopic signa- 
tures can be avoided. The 09dc-int model of lrlachinger et al.l 
(J2012T ). as a first attempt in this direction, has too little flux 
during the tail phase when the interaction should have come 
to an end. 

The 09dc-tail model has been constructed so as to 
improve on this aspect and reproduce the luminosity of 
SN 2009dc during the tail phase between 60 and 180 d af- 
ter the explosion. It has been built with a physical pic- 
ture in mind that is quite similar to that of 09dc-int: the 
kinetic energy of 1.2 x 10 51 erg is typical of a normal 
SN la from a Chandrasekhar-mass WD progenitor. Also 
the 56 Ni mass of 1 .0 Mq is within the ra nge normally ob- 
served in SNe la (|Stritzinger et alj|200q ). though on the 
upper end. While the 09dc-int model only represents the 
ejecta of a Mch WD, the total ejecta mass of 09dc-tail 
is significantly larger. The additional mass can be inter- 
preted as swept-up circumstellar material, leading to re- 
duced ejecta velocities and higher densities compared to a 
freely expanding SN la. At the same time, a small frac- 
tion of the kinetic energy would be transformed into light in 
the ejecta-CSM interaction, potentially compensating the 
flux deficit of the 09dc-tail model at early phases. Whether 
or not this may work and whether the resulting early-time 
spectra would resemble those of superluminous SNe la can 
only be deci ded in detailed radiation-hydrodynamics sim- 
ulations (e.g. iFrver et al.ll20irj ; iBlinnikov fc Sorokinalkuld : 
iNoebauer et al.l I2012T ). which go beyond the scope of this 
work. 

Nevertheless, a few estimates on the required CSM 
properties and a possible progenitor system can be made. 
The CSM would have to be H- and probably He-free, since 
these elements are not detected in the spectra of superlu- 
minous SNe la at any phase. Most likely it would instead 
consist of a mix of C and O, which would be naturally ex- 
plained if it had been produced in a merger of two C/O 
WDs, with the secondary being disrupted and its material 
being accreted slowly on to the primary. The masses of these 
two WDs would have to sum up to ~ 2 Mq . At the time the 
accreting primary approaches Mch and explodes, it would 
be surrounded by 0.6-0.7 Mq of CSM, distributed in an 
extended spherical envelope support ed by thermal pressure 
jShen et al.ll2012l ; ISchwab et al.ll2012l)F|Thi s scenario is very 
similar to that studied bv lFrver et al.l ( 2010T 1. who found that 
in a SN la enshrouded by a C/O-rich CSM originating from 
a merger event the light curve and spectra are altered sub- 



4 We emphasise that our 09dc-tail model has been constructed 
in a self-consistent way, in the sense that the kinetic energy of 
the ejecta agrees with the energy release by nuclear burning (see 
Table [4] for the composition) minus the binding energy of a Mch 
WD. 
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Table 4. Density profile and clement abundances of 09dc-tail. The total mass of the model, 2.06 Mq, is distributed to C (0.24 Mq), 
O (O.47M ), Si (0.25 M Q ), S (0.05 M ), Ca (0.01 M ), stable Fe (0.05 M ) and 56 Ni (0.99 Mq). 



Shell 


"out 


log(p) 6 


M c 

cum 


El 1 

k,cum 


X(C) 


X(O) 


X(Si) 


X(S) 


X(Ca) 


X(stable Fe) 


X( 56 Ni) 


1 


1000 


0.72 


0.011 


6.6 x 10" 5 


0.000 


0.000 


0.000 


0.000 


0.000 


0.050 


0.950 


2 


2000 


0.68 


0.082 


1.9 x 10" 3 


0.000 


0.000 


0.000 


0.000 


0.001 


0.049 


0.950 


3 


3000 


0.64 


0.256 


0.014 


0.000 


0.000 


0.000 


0.000 


0.002 


0.048 


0.950 


1 


4000 


0.52 


0.514 


0.046 


0.000 


0.000 


0.000 


0.000 


0.005 


0.045 


0.950 


5 


5000 


0.32 


0.783 


0.101 


0.000 


0.000 


0.000 


0.000 


0.005 


0.045 


0.950 


6 


6000 


0.12 


1.035 


0.178 


0.020 


0.100 


0.100 


0.020 


0.005 


0.035 


0.720 


7 


8000 


-0.18 


1.447 


0.386 


0.100 


0.260 


0.400 


0.080 


0.005 


0.005 


0.150 


8 


10000 


-0.58 


1.718 


0.608 


0.236 


0.500 


0.200 


0.040 


0.002 


0.002 


0.020 


9 


12000 


-0.98 


1.878 


0.804 


0.318 


0.620 


0.050 


0.010 


0.001 


0.001 


0.000 


10 


14000 


-1.38 


1.967 


0.955 


0.393 


0.600 


0.005 


0.001 


0.000 


0.001 


0.000 


11 


16000 


-1.78 


2.015 


1.062 


0.499 


0.500 


0.000 


0.000 


0.000 


0.001 


0.000 


12 


18000 


-2.18 


2.039 


1.132 


0.499 


0.500 


0.000 


0.000 


0.000 


0.001 


0.000 


13 


20000 


-2.58 


2.051 


1.175 


0.499 


0.500 


0.000 


0.000 


0.000 


0.001 


0.000 


11 


24000 


-3.18 


2.059 


1.214 


0.499 


0.500 


0.000 


0.000 


0.000 


0.001 


0.000 



a Outer boundary of the shell (kins 1 ). b p in gem 



Cumulative mass in Mq. d Cumulative kinetic energy in 10 51 erg. 



stantially by shock heating. The light-curve peak tends to 
become much broader and a strong UV flux is generated, 
a feature that is also o bserved in early spectra of superlu- 
minous SNe la (see e.g. IScalzo et al.ll201(J : I Silverman et al.l 
l201ll ; lTaubenberger et al.ll201ll ). 



To reconcile the synthetic bolometric light curve of 
the 09dc-tail model with observations, additional luminosity 
from interaction is required during the first ~ 60 d after the 
explosion (Fig. [9l bottom panel). This coincides with the 
duration of the peak phase, which may suggest that there is 
not ongoing CSM interaction over this full period of time, 
but that the photons are produced during an early and com- 
paratively short interaction episode, trapped in the optically 
thic k eject a and release d on photon-diffusion time scales (see 
also IScalzo et all l2012h . 

Having swept up 0.6-0.7 Mq of material composed of 
up to 50 per cent of carbon, the ejecta of such an explo- 
sion would contain at least an order of magnitude more car- 
bon than ordinary dela yed detonations of a naked Mch WD 
(|Seitenzahl et alj|2013l ). The strong and persistent Cn fea- 
tures in the early spectra of superluminous SNe la would 
in this scenario not come as a surprise. Moreover, the same 
carbon might form dust at later phases, as already discussed 
in Section 15.21 and thus explain the luminosity drop in at 
least some superluminous SNe la during the nebular phase. 

Of course, the proposed scenario is not free of problems. 
For example, it is not clear how much fine-tuning is neces- 
sary to produce a CSM with the right pro perties in a WD 
merger. As mentioned bv I Fryer et al.l (|2010l ). besides the to- 
tal mass also the density profile of the circumstellar ma- 
terial has an enormous impact on the resulting light curves 
and spectra. To explain observed superluminous SNe la, one 
would need a configuration that boosts the light-curve peak 
for ~ 60 d, slows down the light-curve evolution, but modifies 
the spectra only moderately, preserving most characteristic 
SN la features. Another critical point is that slow mergers of 
such massive WDs might not exist. The closer the mass ratio 
of the primary to the secondary WD is to unity, and the more 
massive the two WDs are, the more likely the merger will 



proceed violently (jPakmor et al.ll2012j ). not leading to an 
extended C/O envelope. Binary population synthesis calcu- 
lations do not pred ict primary WD m asses above 1.3 Mq in 
C/O- WD mergers (|Ruiter et al.ll2013l ). To arrive at a system 
mass of ~ 2 Mq , both WDs hence need to be relatively mas- 
sive. Finally, if non-violent mergers of C/O WDs do exist, 
the primary may collapse to a neutron star owing to electron 
captures rather than explode as a SN la (|Saio fe Nomotd 
ll985l ; IShen et al.ll2012T ). Nonetheless, the proposed scenario 
of a slow merger of two massive C/O WDs that leads to a 
Mch SN la explosion enshrouded by a C/O-rich CSM re- 
mains a promising explanation for superluminous SNe la. 
Its biggest strength is its ability to describe many thus far 
seemingly decoupled properties of superluminous SNe la in 
a unified way as logical consequence of an enshrouded explo- 
sion, with no need to resort to unreasonable 56 Ni and ejecta 
masses. 



6 CONCLUSIONS 

The late-phase photometry and spectra of a group of 
proposed super-Chandrasekhar SNe analysed in this work 
have highlighted interesting trends. The superluminous 
SNe 2006gz, 2007if and 2009dc are distinguished from nor- 
mal and 91T-like SNe la by nebular spectra with very weak 
[Fein] lines and a likely contribution of [Call] to the emis- 
sion around 7300 A. The ejecta of these objects are appar- 
ently not as highly ionised as in normal SNe la, which could 
be a consequence of higher densities due to rather low ex- 
pansion velocities. The proposed super-Chandrasekhar ob- 
ject SNF20080723-012 does not share these characteristics, 
and is instead very similar to SN 1991T in both its nebular 
spectra and its light-curve evolution. We are hence inclined 
to consider SNF20080723-012 as a classical 91T-like object 
rather than postulate fundamental diversity in the late-time 
behaviour of superluminous SNe la. In fact, some diversity 
is present also in the late-time spectral energy distribution 
and light-curve decline of SNe 2006gz, 2007if and 2009dc, 
with two of these objects showing a light-curve break after 
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~ 150-200 d leading to a more rapid lading thereafter. This, 
however, may be understood in terms ol different time scales 
and intensities ol dust formation. 

Accepting this explanation, studying the light-curve tail 
before the possible break (i.e., at phases when the light 
curves are presumably powered by radioactive decay ol 56 Co 
and opacities to optical photons are low) may be a promis- 
ing way to constrain 56 Ni and total ejecta masses ol super- 
luminous SNe la, and thus get a handle on the nature ol 
these objects. To this end we have compared the observed 
UBVRIJ HK-bolometric light curve ol SN 2009dc with syn- 
thetic bolometric light curves ol models proposed in the liter- 
ature. We find that models that have enough 56 Ni to explain 
the light-curve peak by radioactive decay and at the same 
time enough mass to avoid too high ejecta velocities, are al- 
most inevitably too luminous on the radioactive tail because 
ol too strong 7-ray trapping. This dilemma may be overcome 
with models that assume additional luminosity from ejecta- 
CSM interaction during the peak phase. Improving on earlier 
work, we have presented one such toy model with ~ 1 Mq ol 
56 Ni and ~ 2 Mq ol ejecta that provides a convincing match 
to the observed light-curve tail ol SN 2009dc. This model 
may be interpreted as a 'normal' SN la explosion ol a Mch 
WD enshrouded by 0.6-0.7 Mq ol C/O-rich material that is 
swept up as the ejecta expand. Such a configuration could 
be the outcome ol a non-violent merger ol two massive C/O 
WDs. 
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Table Al. Magnitudes of the local sequence stars in the field of 
SNF20080723-012 (Fig.Q] top). 



ID 


B 


V 


R 


I 


1 


20.97 ± 0.03 


20.49 ± 0.02 


20.16 ±0.02 


19.84 ±0.02 


2 


20.64 ±0.03 


19.15 ±0.01 






3 


21.53 ±0.02 


21.01 ±0.03 


20.68 ±0.03 


20.34 ±0.02 


1 


22.05 ±0.02 


21.55 ±0.02 


21.24 ±0.02 


20.86 ±0.03 


5 


21.55 ±0.03 


20.54 ±0.02 


19.87 ±0.03 


19.36 ±0.05 


6 


19.26 ±0.04 


18.81 ±0.04 






7 


19.52 ±0.02 


18.92 ±0.02 






8 


23.06 ±0.03 


21.56 ±0.03 


20.40 ±0.02 


19.03 ±0.06 


9 


22.66 ±0.03 


21.72 ±0.02 


21.10 ±0.04 


20.55 ±0.03 


10 


22.59 ±0.04 


22.09 ±0.02 


21.80 ±0.05 


21.42 ±0.06 


11 


22.00 ±0.04 


20.56 ±0.02 


19.57 ±0.02 




12 


20.32 ±0.04 


19.28 ±0.01 
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Table A2. S- and if -corrections applied to the photometry of SN 2007if and SNF20080723-012. 



MJD 



SN 



Epoch" 



S B 



S v 



Sr 



Si 



K b 



Ki 



Kf. 



Kj Telescope 



54 733.2 


SN 2007if 


358.2 




-0.038 








0.322 






VLT 


54 739.8 


SNF 


54.8 


0.042 


-0.051 


-0.004 




-0.261 


-0.161 


0.102 




CAHA 


54 882.2 


SNF 


187.3 


0.051 


0.016 


0.009 




-0.481 


0.273 


0.223 




CAHA 


54 947.4 


SNF 


248.0 


-0.087 


-0.086 


-0.198 


-0.004 


-0.592 


0.490 


0.283 


0.058 


VLT 


54 973.2 


SNF 


272.0 


-0.087 


-0.086 


-0.198 


-0.004 


-0.592 


0.490 


0.283 


0.058 


VLT 



a Phase in rest-frame days with respect to S-band maximum [MJD 
for SNF20080723-012 jScalzo et alj|2012l) 1. 



54 348.4 for SN 2007if l|Scalzo et alj|2010h and MJD = 54 680.9 



